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ABSTRACT 

We present the preliminary results of a new survey of radio sources using the Ryle 
Telescope at 15.2 GHz. This is the highest frequency at which a survey has been done 
that is relevant to the issue of radio source contamination in CMB experiments. The 
differential source count of the 66 sources found in 63 deg 2 is 80(S/Jy)~ 2,0 Jy~ 1 sr~ 1 , 
from « 20 mjy to « 500 mJy. Extrapolating this to 34 GHz (where many cm-wave 
CMB experiments operate) gives an estimated temperature contribution of sources 
AT con f = 9/iK in a CMB image, with a beam corresponding to multipole I » 500. A 
means of source subtraction is evidently necessary, otherwise the signal-to-noise ratio 
in CMB images will be limited to 4 or 5, becoming worse at higher resolution. We 
compare the population of sources observed in this new survey to that predicted by 
extrapolation from lower frequency surveys, finding that source fluxes, and indeed the 
existence of many sources, cannot be determined by extrapolation. 
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1 INTRODUCTION 

A fundamental limitation to all measurements of primordial 
Cosmic Microwave Background (CMB) anisotropies is that 
set by astronomical foregrounds, both from the Galaxy and 
from extragalactic sources. The degree to which extragalac- 
tic radio sources contaminate is dependent on both the res- 
olution and observing frequency of a particular experiment. 
For example, owing to relatively low collectin g area and an- 
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Since sources are distributed randomly on the sky, the power 
spectrum of point sources also increases as I 2 . The correct 
removal of point sources therefore becomes more important 
as the CMB power spectrum is investigated on smaller an- 
gular scales, a regime in whic h there is even greater int erest 
since the re cent Boomerang ( le Bernardis et al. 2000) and 



gular re solution, the COBE / DM R results ( a.g. Smoot et 
al. 1992 ) are essentially unaffected by extragalactic sources. 
Yet as the new generation of CMB experiments come on- 
line and strive to reach higher resolution and sensitivity, 
foreground sources are becoming an increasingly important 
consideration. 

This can be seen simply by considering the Rayleigh- 
Jeans approximation for the flux density of an unresolved 
source as a function of the angular resolution, D." 1 ^ 2 , of an 
experiment: 



MAXIMA flHanany et al. 2000D results. 

The removal of foreground sources, however, is not triv- 
ial. They are known to have a wide range of spectral indices. 
At low frequencies, their spectral indices a (with 5*oc v~ a ) 
are usually steep (a > 0.5), but some are flat (a = 0) or 
may even be rising. Very occasionally, a spectrum may rise 
s o steeply that it mimics that of the CMB with a — — 2 
(lEdge et al. 19981). 



S = 



2k B TD. 
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(1) 



where S is the source flux density, fee is the Boltzmann con- 



The secure method of determining the point sources 
present in a CMB observation is to map at high resolu- 
tion. CMB experiments, however, operate at frequencies (> 
30 GHz) where the population of point sources is not well 
characterised. Indeed, until the present results, ther e existed 
no suitable high -frequency survey a bove 4.8 GHz (Gregory 



stant, f ! is the beam solid angle, A the observing wavelength & Condon 1991). Windhorst et al. (1993) provide counts at 



and T is the resulting antenna temperature. Expressing this 
in terms of multipole index / (where I ~ 2n/0 for high I) we 
see that the temperature contribution of a source increases 
as I 2 : 



8.4 GHz of the sub-mjy source population, but these are too 
faint to affect primordial C MB stu dies. They also plot the 
10.0-GHz data of Aizu et al. ( |l987| ), spanning 30-3000 mjy. 
However, as Aizu et al. state, their 10.0-GHz observations 
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only include sources selected at 5 GHz and do not represent 
an independent survey. Theoretical predictions of the extra- 
galactic source c ounts at a wide rang e of frequencies have 
also been made (Foffolatti et al. 1998), but observations at 
relevant frequencies are still needed to confirm the results. 
Furthermore, although source counts enable one to predict 
the point source contribution to a measured CMB power 
spectrum, this is helpful only if the source contribution is 
small compared with the CMB power. In general, finding 
the sources with high-resolution imaging is necessary. 

In this paper we present the preliminary results of 
a source survey at 15 GHz that is part of the source- 
subtraction strategy for a new CMB interferometer, the Very 
Small Array (VSA, see Jones & Scott 1998). This is the 
highest frequency at which a survey has been undertaken 
that is relevant to the issue of radio-source contamination in 
cm-wavelength CMB experiments. The spectra and number 
count of the sources are discussed, as are the implications 
for the removal of extragalactic sources from CMB measure- 
ments. We stress that the temperature confusion estimates 
presented here apply to images of the CMB; we shall present 
estimates of the effects of sources on the CMB power spec- 
trum in a subsequent paper. 
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Figure 1. The differential source count of sources observed in 
VSA1 at 15.2 GHz. The fall-off at low S is ascribed to the lack of 
completeness. 



taneously with VSA observations so that variable sources 
can be subtracted accurately. 



2 SOURCE SUBTRACTION FOR THE VSA 
AND THE FIRST 15-GHZ SURVEY 

The VSA is a close-packed, 14-element interferometer for 
sensitive measurement of CMB anisotropies. In its compact 
configuration the VSA will measure the CMB power spec- 
trum in the Z-range 120-700, whilst in its extended array it 
will reach I — 1800. The telescope is currently operating at 
34.1 GHz and we expect extragalactic radio sources to be a 
major contaminant in VSA observations. We have developed 
a strategy to deal with this problem in two steps. 

First, we use low - frequency surveys, NVSS at 1.4 GHz 
( Condon et al. 1998 ) and Green Bank at 4.8 GHz ( |Grcgory 



fa Condon 1991), to select CMB fields in which there are 
predicted to be no sources brighter than 500 mjy at 34 GHz. 
Predictions are made by extrapolating the fluxes of every 
source in the 4.8-GHz catalogue to 34 GHz on the basis 
of its spectrum between 1.4 and 4.8 GHz. The VSA fields 
are also selected to be as free as possible of other Galactic 
and extragalactic foregrounds. In particular they have been 
chosen for high galactic latitude, low galactic synchrotron 
and free-free emission, low dust and putative dust, and a 
lack of nearby large-scale structure. 

Having selected suitable fields, and prior to VSA obser- 
vation, we locate potentially contaminating radio sources 
by mapping each VSA field at a much higher resolution 
than that of the VSA. This is achieved using the Ryle Tele- 
scope (RT) in C ambridge, operating at 15 .2 GHz in a raster- 
scanning mode ( Waldram & Pooley 2000 ) . We reach an rms 
noise level of a — 4 mjy in order to identify all sources 
above 20 mjy at 15.2 GHz. Simultaneously with VSA ob- 
serving, we then monitor the flux of each source using a sepa- 
rate single-baseline interferometer operating at the same fre- 
quency as the VSA. This interferometer uses two 3.7-m an- 
tennas located next to the VSA, using identical receivers and 
correlator to the main array. The monitoring is done simul- 



3 RESULTS AND IMPLICATIONS OF THE 
15-GHZ SURVEY 

The first VSA field, VSA1, covering an area of 63 deg 2 and 
centred at 00 ,I 20 m -|- 30° (B1950), contains 66 sources with 
S <; 20 mjy. The differential source count for this field is 
shown in Figure 1 and, between 20 and 500 mjy, is well 
fitted by 



n(S) = 8 
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The spectral index distribution between 1.4 GHz and 
15 GHz for all the sources found in our survey is shown 
in Figure 2. A significant fraction of these sources were not 
predicted to be seen at all; they have components with flat 
or rising spectra that only surface above 4.8 GHz. Note that 
in order to assign a spectral index to all the sources, we have 
had to use a\ 5 4 rather than a\% since the Green Bank 4.8- 
GHz survey is far less deep than NVSS. A consequence of 
this is that a\ 5 8 is probably even steeper than a\ 4 for a few 
of the sources. 

To estimate th e effe ct of these sources on CMB images, 
we follow Scheuer (1957): 



n(S)S 2 dS, 



(4) 



where er con f is the confusion noise in flux density from unsub- 
tracted sources and Sn m is the flux density above which all 
sources have been subtracted. For I —500 (typical of the VSA 
compact array and of other cm-wave experiments) and as- 
suming a Gaussian distribution of baselines, S7 = (27r/500) 2 
giving a 43-arcmin beam. If no source subtraction is at- 
tempted (i.e. Sum m 500 mjy, and we emphasize here that 
this VSA field has been pre-selected to avoid bright sources) , 
then equation (4) predicts a con f — 79 mjy at 15 GHz. 
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Figure 2. Spectral index distribution, ctf 5 4 , of sources found at 
15.2 GHz with the Ryle Telescope. Grey bars represent observed 
sources predicted by extrapolation of spectra using the 1.4-GHz 
NVSS and 4.8-GHz Green Bank surveys. The black bars represent 
those sources observed but not predicted. 



We now consider what cr con f might be at the VSA ob- 
serving frequency of 34 GHz. Our educated guess of the 
mean spectral index between 15 and 34 GHz, given the 
change in source population from 4.8 to 15 GHz that we 
have found at these flux levels, is that af| = 0.5. This is less 
steep than the typical a at low frequencies but reflects the 
increasing fraction of flat- and rising-spectrum sources at 
high frequency. Assuming this value of a, we estimate CT con f 
at 34 GHz to be 52 mjy in a 43-arcmin beam correspond- 
ing to I S3 500. Equation (1) then implies that the resulting 
temperature from unsubtracted sources, AT con f, is 9.2 fjK. 
This would limit the signal-to-noise ratio in CMB images 
to 4 or 5, and the problem becomes even worse at higher 
resolution. Source subtraction is evidently necessary. 

By surveying for sources one can obviously decrease 
Sii m . If a survey reaches Su m = 50 mjy at 15 GHz, then 
Equations (1) and (4) and the assumed af| = 0.5 give 
AT con f = 2.9 /uK at 34 GHz in a beam corresponding to 
I S3 500. The fact that the RT survey reaches S3 20 mjy 
means that it will find all sources which have af| rising as 
steeply as — 1. This is likely to include the majority of the 
source population. Sources with a approaching —2, however, 
are known to exist (see e.g. Edge et al. 1998). An upper 
limit to the contribution of such sources can be found as fol- 
lows. A source with af| = —2 and S15 = 20 mjy will have 
S34 = 103 mjy. If we assume that all sources have af| = —2 
and that n(S) is given by Equation (3), then the RT survey 
effectively has Sn m = 103 mjy at 34 GHz. This results in 
AT con f = 6.4 /iK in a beam corresponding to I S3 500. The 
population will of course not be this extreme, and AT con f 
for VSA compact array maps is likely to be close to 3 /iK. 

Finally, we consider the effect of source subtraction us- 
ing only predicted source fluxes and positions extrapolated 
from low frequency surveys. Using this approach we expect 
to find 122 sources with a flux density greater than 20 mjy 
within our survey region of 63 deg . We find, however, that 
only 55 of the predicted sources are actually found, the rest 
having spectra that fall off enough between 4.8 and 15 GHz 
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Figure 3. The differential count of residual flux remaining after 
source subtraction has been done using predicted source fluxes. 



to drop them below our flux limit. Furthermore, 11 sources 
were found in our survey which were not predicted by ex- 
trapolation. Five of the unpredicted sources do not even 
appear in the 4.8-GHz Green Bank catalogue. To find the 
residual source flux that would be left on an image if only 
predicted source fluxes had been subtracted, we calculate 
the modulus of the difference between predicted and ob- 
served flux for each source. Figure 3 shows the differential 
count of this residual flux; note that the remaining source 
contamination is somewhat worse than if no source subtrac- 
tion had been attempted. Specifically, we find that substan- 
tial errors in individual source fluxes are made, both in over- 
predicting the fluxes of a large fraction of sources with down- 
curving spectra, and non-prediction of a significant fraction 
of sources with rising spectra. This result, however, is not 
unexpected. As we move to higher frequency regimes, we 
expect many more sources to be synchrotron-self-absorbed 
(SSA), exhibiting intrinsically rising spectra and increased 
variability. Since the surveys we used to predict the source 
population at 15 GHz are a decade lower in frequency, one 
would therefore expect to underestimate the number of SSA 
sources, leading to a non-prediction of many sources with 
rising spectra. Furthermore, the surveys used were under- 
taken at different epochs and hence no allowance can be 
made for variable sources. 



4 CONCLUSIONS 

(1) We have used the Ryle Telescope at 15.2 GHz to sur- 
vey 63 deg 2 of sky chosen, on the basis of low-frequency in- 
formation, to contain no source with flux density > 500 mjy. 
The differential source count of the 66 sources found is 



80(Syjy)- 2 -°Jy- 



from « 20 mjy to S3 500 mjy. 



(2) Extrapolating the source count to 34 GHz, near to 
which current cm-wave CMB experiments operate, the tem- 
perature contribution of sources to a CMB image with a 
43-arcmin beam (corresponding to I S3 500) is estimated to 
be ATconf = 9/iK. Hence, the signal-to-noise ratio in CMB 
images will be limited to 4 or 5, and even more at higher 
resolution, if foreground sources are not removed. 
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(3) In the case of the VS A compact array, which will mea- 
sure the CMB power spectrum in the /-range 120-700, the 
RT source surveying should reduce AT con f to 3 /iK in a 43- 
arcmin beam at 34 GHz. 

(4) Our results show that the existence and flux density 
of individual sources cannot be predicted accurately by ex- 
trapolation from low-frequency surveys. We find that trying 
to predict source fluxes at 15 GHz by extrapolating their 
1.4 and 4.8-GHz fluxes gives a greater error, i.e. a greater 
ATconf, than attempting no source subraction at all. 
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